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1. Introduction 

The Atacama Large Millimeter/submillimeter Array (ALMA)* is the largest telescope yet con¬ 
structed. Initial results from its early science observations span the range from the detection of a 
remarkable spiral patterned mass loss resulting from the interaction of an ABG star and its binary 
companion (Maercker et al. 2012) to constraining the density of extreme starburst galaxies at z > 1 
(Karim et al. 2013). Here we focus on galactic star formation and the synergies provided by the 
complementarity of ALMA and SKA observations. 

For studies of star formation, ALMA observations will trace how cold dust and gas is assem¬ 
bled on the small-scales within a molecular cloud and how disks around young stars are formed. 
SKA will be essential for complementing these studies by tracing the higher energy phenomena 
associated with the accretion of this material on to the star, probing the stellar and disk magne¬ 
tospheres, and the ionising feedback from outflows and stellar winds. Combined, these cutting 
edge facilities will reveal the process of mass assembly and accretion onto young stars and will be 
revolutionary for studies of star formation. 


2. Capabilities: ALMA and SKA 

Comprising fifty four 12m diameter dishes plus twelve 7m diameter dishes, ALMA is an 
international collaboration between ESC, NRAO and NAOJ presenting Europe, North America 
and East Asia respectively in cooperation with the Republic of Chile. Operating between 30 GHz 
and 900 GHz (10 mm - 0.35 mm) in 10 frequency bands which span the available atmospheric 
windows to study the emission from cold dust and molecular gas, AEMA is located on the 5000m 
altitude Chajnantor plateau in northern Chile. 

There were two key science drivers for the design of AEMA. The first is the imaging of the 
physical, chemical and magnetic structure of protostellar/protoplanetary disks around low mass 
protostars in the nearest star forming regions. The second is imaging of Milky Way-mass galaxies 
at z = 3 in CO or C+ (the brightest tracers of molecular gas) in less than 24 hours. In parallel with 
these is the requirement to routinely deliver high precision images at an angular resolution of 0.1”. 
With 8 GHz of instantaneous bandwidth, the AEMA correlator provides spectral resolutions up to 
maximum of 3.8 kHz, corresponding to a velocity resolution of 0.01 km/s at 110 GHz. In full oper¬ 
ation AEMA with reach continuum noise levels of 1 mJy (at 110 GHz) in 1 second, corresponding 
to a mass of ~0.004 Mq per beam in the nearest star forming regions and ~ 5 Mq at 5 kpc. Eigure 1 
shows the angular size scale versus frequency coverage of AEMA compared with SKAl-Mid (and 
the NRAO Karl G. Jansky Very Earge Array, JVEA). As the figure demonstrates, to probe similar 
angular scales as AEMA at frequencies above 60 GHz, requires frequency Band 5 on SKAl-Mid. 
As an indication of the range of physical scales relevant for studies of galactic star formation, the 
horizontal lines on the figure show the angular size required to image a hypercompact HII region, 
the earliest phase of massive star formation where ionised gas is detected, and a large circumstellar 
disk/toroid at the typical distance to molecular clouds in the molecular ring, ~ 5 kpc. 

* www.almaobservatory.org 
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Frequency (GHz) 

Figure 1: Angular size scale as a function of frequency for SKAl, JVLA and ALMA. Frequency band 5 on 
SKA 1-Mid is identified separately. The shaded regions show the angular scale from the maximum resolution 
at each frequency to the instantaneous field of view. Larger angular size scales are of course accessible 
through mosaicing. The two horizontal lines indicate the resolution required to resolve a 0.01 pc-sized 
hypercompact HII region (upper) and a 1000 AU circumstellar disk/toroid (lower) at a distance of 5 kpc. 


3. Tracing global collapse, disks and outflows 

The ultimate fate of infalling material, and how stars actually gain mass, depends on the 
processes taking place on small spatial scales close to a forming star. Recent evidence suggests 
that the final stages of accretion on to the central star may be highly episodic. For example, 
ALMA observations of the chemical composition of the envelope around a low mass protostar 
shows evidence of heating associated with a burst of accretion 100-1000 years ago when the mass 
accretion rate was a factor of ~ 100 above its current rate (Jprgensen et al. 2013). 

On the large spatial scales, ALMA observations of the massive infrared dark cloud SDC335 
(Peretto et al. 2013) show that the infalling dense gas seen in single dish images is confined fo 
a number of filamenfs, allowing us fo frace fhe large-scale global collapse. On fhe small scales, 
ALMA observafions of high excifafion lines of bofh complex organic species (Pineda el al. 2012) 
and simple species (Zapala el al. 2013) show infall profiles fowards one componenl of fhe low 
mass proloslellar binary IRAS 16293. Indeed, measuring infall is a crucially imporlanf fool for 
underslanding slar formation: line profiles diagnostic of infalling gas (e.g. Anglada el al. 1987) 
have been identified towards slarless cores (e.g. Lee el al. 1999), low mass proloslars (e.g. Walsh 
el al. 2006), high mass young slellar objecls (e.g. Fuller el al. 2005), UCHII regions (e.g. Ho & 
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Young 1996), and molecular filaments (e.g. Kirk et al. 2013). 

Circumstellar disks are the mass reservoir out of which planets form. ALMA observations are 
elucidating the structure and composition of these disks with unprecedented detail. Observations 
of transition disks, relatively evolved disks with larger inner holes, have revealed large azimuthal 
asymmetries in the dust continuum emission at radii of ~ 30 — 50 AU (Casassus et al. 2013; van 
der Marel et al. 2013; Perez et al. 2014). These asymmetries are interpreted as resulting from the 
trapping of dust in large anticyclonic vortices in the disk. 

The growth of dust grains is a critical step towards planet formation which SKA will probe. 
Regions of disks where this growth may be occurring are revealed by ALMA observations, not only 
of the dust continuum emission, but also using molecular lines to trace the CO snow-line, where 
CO condensed on to grain surfaces (Qi et al. 2013; Mathews et al. 2013). ALMA observations 
also reveal evidence for complex chemistry in a transition disk (van der Marel et al. 2014) as well 
the presence of a simple sugars (Jprgensen et al. 2012). Future SKA observations will be critical 
for the identification of complex organic species towards these objects. 

Material in protoplanetary disks can be eroded by the winds and high energy radiation from 
nearby high-mass stars. To quantify this interaction, ALMA observations can probe molecular gas 
in the disk (Mann et al. 2014) while SKA can probe the photoionized gas flow from the disk as well 
as the interaction zone between the gas and the impinging stellar wind (e.g. Graham et al. 2002). 

Outflows are an ubiquitous feature of star formation. Driven by winds and jets which originate 
from the inner disk/disk-star interface region (Li et al. 2014), they not only sculpt, and eventually 
clear the material around the protostar, but also inject energy into the surrounding medium and 
remove angular momentum. The bulk of the molecular line emission seen from outflows detected 
via ALMA observations (e.g. Arce et al. 2013) will trace material swept up by the stellar wind. 
The nature of the jets from stars evolves as the driving source evolves (Frank et al. 2014) but 
observations show the jets have a significant atomic component potentially observable by SKA in 
HI. 

4. Accretion and Variability in Young Stars 

Low-mass young stellar objects (YSOs) are thought to accrete mass via their circumstellar 
disks. Associated accretion shocks can be observed both when mass is accreted onto the disk and 
also when it is eventually accreted onto the central object (e.g. Hartmann 2009). The accretion 
process is observable at wavelengths ranging from the radio to the X-ray regime, including optical 
and infrared wavelengths. The connection of accretion and high-energy processes (e.g. Feigelson 
& Montmerle 1999) detected in both X-ray and radio emission has been highlighted by spectacular 
examples like VI647 Ori (Kastner et al. 2004), while spectral lines in the X-ray can provide an 
estimate of the the accretion rate (Brickhouse et al. 2012). 

In the centimeter regime, accretion processes can be associated with both thermal (free-free) 
and non-thermal (gyrosynchrotron) radiation, as discriminated by their spectral and polarization 
properties. Generally, it is thought that the non-thermal radiation emanates from the magneto- 
spheric structures in the innermost vicinity of YSOs while thermal emission most likely traces the 
bases of outflows and jets further away from the central object. 
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Recently there has been considerable observational progress in two relevant areas. First at cen¬ 
timeter wavelengths, new instrumentation, provided via the expanded capabilities of the JVLA, is 
rekindling protostellar radio astronomy. In tandem with these advances, a multi-wavelength, time- 
domain view of the dynamic processes in YSOs, including accretion is emerging. Well sampled 
time series datasets are now becoming available, in particular in the infrared. Simultaneously ob¬ 
tained multi-wavelength time series datasets are also becoming increasingly important to constrain 
the underlying physical processes. 

4.1 Identifying radio counterparts to low-mass YSOs 

Observations with the VLA initiated protostellar radio astronomy by providing excellent an¬ 
gular resolution to reliably identify radio counterparts to low-mass YSOs. While the first radio- 
detected YSOs identified by the VLA were found in the Orion Nebula Cluster (e.g. Moran et al. 
1983; Garay et al. 1987), the radio sample remained incomplete, as indicated by the fraction of 
radio detected YSOs being consistently lower than the X-ray detection fraction (Forbrich & Wolk 
2013). 

Prior to the sensitivity upgrade of the VLA, only the most nearby star-forming regions could 
be studied in detail. The low-mass Coronet cluster at a distance of 130 pc is an example where 
almost all known YSOs in the inner cluster have radio counterparts (Brown 1987; Forbrich et al. 
2006). This region was shown to host the first of only a few known protostars with confirmed 
non-fhermal radio emission as inferred from circular polarizafion indicative of gyrosynchrofron 
emission (Feigelson el al. 1998; Choi el al. 2009). 

The newly upgraded JVLA has already significanlly advanced Ibis research by providing a 
more complele census of YSOs delected al cenlimelre wavelenglhs (Dzib el al. 2013b; Kounkel 
el al. 2014). Indeed, wilh Ihe improved JVLA, a more complete census of YSOs in nearby slar- 
forming regions is emerging, bolh in terms of sensilivily and in area covered. Using rapid vari- 
abilily, negative speclral indices and polarization as indicators, a recenl, large survey of p Oph 
revealed lhal aboul half of Ihe identified YSOs show non-lhermal emission (Dzib el al. 2013b). 

Similarly, a recenl large-scale (more lhan 2 square degrees) survey of Ihe Orion Nebula Clus¬ 
ter, detected 374 sources al 4.5 GHz (Kounkel el al. 2014). Of Ihese, 148 had been previously 
classified as YSOs and 86 additional sources are inferred to be new YSO candidates. Wilh reliably 
determined speclral indices, Ihese sources will be used as largels for follow-up VLBI observations 
to sludy parallaxes and proper motions. New and sensitive monitoring observations of Ihe Coronet 
cluster show lhal Ihe youngesl sources in Ihe sample are Ihe brighlesl and leasl variable, possibly 
exhibiting moslly Ihermal wind emission (Liu el al. 2014). 

4.2 Time variability 

The signilicanl increase in sensitivity provided by new JVLA observations is also slarling to 
produce radio lighl-curves for YSOs. While YSO X-ray flares have been known for some time (e.g. 
Gelman el al. 2008), very little is known aboul protostellar radio flares (Bower el al. 2003; Forbrich 
el al. 2008) and a possible X-ray-radio connection (e.g. Guedel & Benz 1993; Forbrich & Wolk 
2013). As a resull Ihe physics of YSO radio flares is nol currenlly well understood allhough Ihey 
are Ihoughl to be produced by coronal-lype activity in scenarios lhal also produce X-ray emission 
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(e.g. Drake et al. 1992). Some of the aetivity eould also be direetly or indireetly related to aeeretion. 
However, some YSO flares are thought to be due to the interaetion of large magnetospherie features 
in elose binary systems (e.g., V773 Tau: Massi et al. (2006) and DQ Tau: Salter et al. (2008); 
Getman et al. (2011) and referenees therein) or from periodie aeeretion bursts (Muzerolle et al. 
2013; Balog et al. 2014; Bary & Petersen 2014). 

Impressive time domain studies are now being earried out in the infrared. A 30-day Spitzer 
and CoRot photometrie monitoring eampaign of more than a thousand YSOs in NGC 2264 (Cody 
et al. 2014) elassified a variety of different variability meehanisms in YSOs, mainly on timeseales 
of days. Other reeent work suggests that short-term optieal variability may be due to enhaneed 
aeeretion aetivity (Stauffer et al. 2014), while mid-infrared variability may be due to struetural 
perturbations in the inner disk (Flaherty et al. 2013). Speetroseopie monitoring eampaigns have 
also revealed a new level of eomplexity. Multi-epoeh near-infrared speetroseopy of several ae- 
ereting YSOs in p Oph observed as part of the Spitzer YSOVAR projeet, showed no eorrelation 
between the YSO mid-lR light curves and time-resolved veiling or mass accretion rates (Faesi et al. 
2012 ). 

4.3 Accretion onto high-mass stars 

Strong magnetic fields and circumstellar disks are key drivers of the evolution of, and activity 
in, the inner circumstellar regions of low mass stars and there is growing evidence these impor¬ 
tant components are also present in young high-mass, M > SM©, stars. Keplerian disks have been 
detected through high angular resolution imaging of molecular lines around a number of B-type 
stars (e.g. Cesaroni et al. 2014; Beltran et al. 2014) where, like their lower mass counterparts, they 
presumably play a role in the accretion of material on to the central star. Stellar magnetic fields in 
the range of hundreds to thousands of kG have also been measured towards OB stars (Naze 2014). 
Since stellar magnetic fields are expected to decay away as the high-mass star evolves due to the 
absence of a convective zone, similar or stronger fields are likely present towards younger stars 
where they can play an important role in the rotational braking of the stars, setting their initial rota¬ 
tion rates (Meynet et al. 2011; Rosen et al. 2012). The presence of disks and strong magnetic fields 
suggest that young high-mass stars can show a similar range of circumstellar energetic phenomena 
as seen towards young low-mass stars, probing the final stages of mass inflow on to the central star 
as well as the launching of outflows. 

Indeed non-thermal radio emission has been detected towards an increasing number of young 
high-mass stars (e.g. Andre et al. 1988; Rodriguez et al. 2012; Moscadelli et al. 2013; Rodriguez 
et al. 2014). This emission can arise from a range of phenomena including magnetospherie activity 
(such as due to variations in magnetically mediated accretion on to the central star), interacting 
winds in a massive binary system and synchrotron emission from a jet, all processes which are 
likely to be variable on a range of timescales as seen towards low mass sources (e.g. Dzib et al. 
2013b). For example, a source in Orion has been seen to increase in flux by factor of 3 increase 
over 1 hour (Gomez et al. 2008). 

Thermal radio emission from young high-mass stars arises from their ionised winds and their 
self-photoionised disks, envelopes and surroundings in the form of Hll regions. The accretion 
flow into a Hll region and on to forming massive stars will change the ionization balance and 
hence, flux and size of the ionised region. Such inflowing material will be clumpy and so the mass 
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accretion variable. This is borne out by simulations which show that the accretion rates onto stars 
has considerable time variability, with peak rates between 10 and 100 times higher than the time- 
averaged values (e.g. Klassen et al. 2012b). This variable inflow results in a flickering of the radio 
continuum emission of the HII regions around massive stars. The simulations by Galvan-Madrid 
et al. (2011) found that about 10% of HII regions showed flux variations of 10% or more in 10 
year timescales. A similar result was found by Klassen et al. (2012a). The radio variability of HII 
regions can therefore probe the structure and variability of the accretion into massive stars. 

To date there have only been a very limited number of studies of the radio variability of HII 
regions. Galvan-Madrid et al. (2008) identified a ~ 45% decrease in flux (at 5 GHz) of a hyper¬ 
compact HII (HCHII) region over a 5 year period due to enhanced accretion, while 10% of the 
ultracompact HII (UCHII) regions with Sgr B2 showed significant changes in flux over a timescale 
of 23 years (De Free et al. 2014). The variable central source in W3(OH) has changed flux by a 
factor of 5 on timescales of 9 years (Dzib et al. 2013a) which is interpreted as due to changes in 
the ionised atmosphere of a circumstellar disk possibly due to changes in the accretion through the 
circumstellar disk. 

The masers observed towards many embedded massive stars provide an additional sensitive 
probe of changes in the continuum emission of their exciting sources. For example, towards some 
sources the 6.7 GHz methanol maser emission varies periodically or quasi-periodically (Goedhart 
et al. 2014; Szymczak et al. 2014) which may reflect variations in either the free-free background 
radiation due to changes in the ionization of the circumstellar material or changes in the infrared 
pump radiation. Both of these effects may arise from periodic accretion from a circumstellar disk. 

5. SKAl Young Star Variability Survey 

Both low- and high-mass young stars show centimetre radio continuum emission which is 
variable on a range of timescales from sub-hour to years and decades. This emission traces a range 
of phenomena associated with accretion, magnetospheric activity, and outflow which link the hot 
inner circumstellar regions with the infall and outflow of cooler material traced on larger scales by 
ALMA. SKA time domain continuum surveys of young stars will provide powerful probes of the 
inner circumstellar regions of young stars, opening new windows on the star formation process. 
For example, providing the first comprehensive studies of the properties of the episodic accretion 
in the inner circumstellar regions of forming massive stars. 

In order to confirm the association of radio emission with a particular YSO in crowded clus¬ 
ter environments, as well as isolate the emission close to the central star, such a SKA survey 
will require sub-arcsecond angular resolution. In addition it will require near simultaneous wide 
multi-frequency coverage (Figure 2) to constrain the spectral index of the emission to discriminate 
between thermal and non-thermal emission as well as optically thick and optically thin emission. 
Measurements of the thermal emission in both the optically thick and thin regimes are important to 
distinguish between changes in size of the HII region and its ionization. Full Stokes synthesis will 
be important not only for confirmation of the presence of non-thermal emission, but also to allow 
the separation of thermal and non-thermal components of the emission. 

Figure 2 shows the spectral flux distribution of young HII regions as well as thermal wind 
and synchrotron emission from lower mass stars. To distinguish between the possible emission 
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Figure 2: Spectral flux distribution for typical UCHII region (solid, blue curve) with an emission measure 
of lO^pccm^®, a typical HCHII region (dashed, red curve) with an emission measure of 10^ pccm^®, a 
typical ionised wind (with the flux scaling as V** ®) (green dots) and synchrotron emission ^ (dotted, 

blue curve). The break in the spectrum of the UCHII and HCHII regions marks the transition between 
optically thin and optically thick emission. The wind and synchrotron emission are normalised to the 4.5 
GHz measured flux of young stellar objects detected in Orion (Kounkel et al. 2014) scaled to a distance of 
5 kpc. The vertical dashed lines indicate the survey frequencies while the horizontal dashed lines shows the 
survey 5-(J noise level (5 /iJy). 


mechanisms, observations at three frequencies are required. A proposed set of survey frequencies, 
1.6 GHz, 7 GHz and 14 GHz, is shown on the figure. As well as providing good sampling of the 
continuum spectrum of the sources, these particular frequencies also cover several of the key maser 
transitions (ground state OH masers at 1.6 GHz, Class II methanol masers at 6.7 GHz and 12.2 GHz, 
and the H 2 CO 2(1,1)- 2(1, 2) transition at 14.5 GHz^) which are themselves variable, in some cases 
periodically, and provide an additional, sensitive tracer of changes in the emission from the central 
sources. A summary of the parameters and the required observing time for a sample survey are 
given in Table 1 . 

To sample the range of physical phenomena which give rise to variable radio emission requires 
repeated observations with a wide range of cadences. Analysis of subsets of the observations will 
sample sub-hour variability due to flares arising from accretion or magnetospheric events (e.g. 
Salter et al. 2008; Massi et al. 2006) and rapidly rotating stars (Wolff et al. 2006) while daily to 

^This line is not currently expected to be covered by SKAl although this would be desirable. However, this fre¬ 
quency is expected to be available with the full SKA. 
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Frequency 

(GHz) 

F 

^ src 

(AtJy) 

a 

(AiJy) 

d 

(”) 

Time 

(hours/field) 

Comments 

1.6 

14 

1 

0.5 

5 

Optically thick UCHll regions; synchrotron; 
OH ground state masers 

7 

9 

1 

0.5 

2.5 

Optically thin UCHll, thick HCHll; 
synchrotron-thermal wind cross over; 
6.7 GHz methanol masers 

14 

14 

1 

0.5 

2.5 

Low optical depth HCHll; thermal wind 
emission; 12.2 GHz methanol masers; H 2 CO 

transition 


Table 1: SKAl young star variability survey parameters. The columns show the frequency of observation, 
the target source flux (Fsrc) and rms flux noise levels (o’), the required angular resolution (0) and the integra¬ 
tion time per field. The comments indidate the emission probed and spectral lines which would be covered 
in the observations. 

monthly observations will sample rotating stars and binary systems with longer periods. Longer ob¬ 
servation intervals will sample processes ranging from outflow bursts to changes in the ionization, 
and therefore flickering, of HII regions due to accretion flows (Galvan-Madrid et al. 2011). 

The angular resolution of the survey is set by the requirement to isolate individual sources in 
crowded star forming regions and resolve the emission from individual HCHll regions with sizes 
~ 0.02 pc out to at least the molecular ring at ~ 5 kpc. To connect the detailed studies which can 
be carried out in nearby star forming regions with the processes in the wider galactic plane requires 
that similar sources are detectable in both kinds of regions. Therefore the sensitivity required for 
each single observation is set by the need to sample a similar range of sources to those observed 
in nearby star forming regions such as Orion at the typical distance of a star forming region in the 
molecular ring, 5 kpc. To detect wind emission sources and synchrotron sources in the Orion star 
forming region (Kounkel et al. 2014) at 5 kpc implies a required noise level of 1 /r Jy at each of the 
observing frequencies (Fig. 2) which should provide detections of up to few hundred sources per 
star forming region. 

Combining the observations at the different epochs will also provide a highly sensitive sur¬ 
vey for faint sources including Hll regions around relatively low mass stars. The combination of 
14 GHz observations at 10 epochs would reach a noise level of ~ 0.3 /f Jy, sensitive at a 6-a level 
to the Hll regions from stars with masses as low as h.SM© corresponding to a spectral type of B5 
(Diaz-Miller et al. 1998). 

6. Summary 

A SKA time-domain survey of star forming regions will produce well sampled radio light- 
curves of hundreds of sources per star forming region. As is the case for infrared light-curves of 
young sources, these will undoubtably have a range of properties, tracing and constraining a range 
of phenomena. For example, the identification of rotationally modulated emission will provide 
the first observational constraints on the rotation rates of embedded protostars. Observations of the 
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flickering of HII regions will trace the inflow of material around the highest mass objects while flare 
activity in these, and other, sources can trace the final infall of material on to the central protostars. 
These observations can not only constrain the properties of the accretion but also the stochastic 
heating of the circumstellar regions which can surpress fragmentation (e.g. Krumholz et al. 2007), 
leading to the formation of more massive stars. 

Looking forward to full SKA the factor 10 enhancement in sensitivity will enable studies of 
smaller (and therefore also shorter time scale) flux variations. Full SKA will also make stellar 
clusters to beyond the Galactic centre accessible to detailed time-domain studies, probing star 
formation over a much wider range of environments. The enhanced resolution of SKA will allow 
studies of the changes in morphology of the sources during their variation in flux, provide stronger 
constraints on the mechanisms involved and the models for the emission. 

Building a comprehensive model for the formation of both low-mass and high-mass stars re¬ 
quires understanding the evolution of gas and dust from molecular clouds down through clumps 
and cores and eventually on to the forming stars. This is only possible with the combination of 
ALMA to trace the cool molecular gas and dust and SKA to follow this material into the inner 
circumstellar regions. The radio light-curves which SKA will produce for hundreds of sources in 
star forming regions out to 5 kpc will provide the first detailed insight across a wide range of stel¬ 
lar masses of the transient energetic phenomena occurring on the small spatial scales close to the 
central star. These light curves will for the first time provide a systematic survey which can study 
the magnetospheric interactions in young binary systems, tracing episodic accretion, the origin of 
outflows and potentially constrain the rotation rates of deeply embedded sources. 
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